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Abstract 

The new, model-independent method to estimate simultaneously the level densities excited 
in the (n, 7) reaction and the radiative strength functions of dipole transitions is developed. The 
method can be applied for any nucleus and reaction followed by cascade 7-emission. It is just 
necessary to measure the intensities of two-step 7-cascades depopulating one or several high- 
excited states and determine the quanta ordering in the main portion of the observed cascades. 
The method provides a sufficiently narrow interval of most probable densities of levels with 
given J* and radiative strength functions of dipole transitions populating them. 



1 Introduction 

The observed parameters of the cascade 7-decay of the compound nucleus can be reproduced in 
the calculation if one determines (in the frameworks of some model) at least 

1. the mean density p of the excited states with given spin and parity J n , and 

2. the mean width T\i of 7-transitions between the arbitrary states A and i. 

The objects of primary interest are the total radiative width T 7 of the compound nucleus (neutron 
resonance) and the spectrum of 7-emission. It may be, for example, the intensity i" 77 of the cascades 
of two successive 7-transitions between the compound state and given low-lying level via a great 
number of intermediate levels. The experimental data on J~~ (as a function of the energy of their 
intermediate level) are obtained for over 30 nuclei from the mass region 114 < A < 200 (see, 
e.g., HD w hh a precision of approximately 10%. The experimental values of T 7 are known within 
the same accuracy. Unfortunately, such accuracy cannot be achieved in the calculation of these 
parameters for an arbitrary nucleus because there are no models that would predict p and T\i with 
the mentioned above precision. 

This is seen from the data of Table 1 which represent the mean ratio between the experimental 
cascade intensities I 77 p and those calculated I 77 using the known enough models of level density 
[§, ||] and radiative strength functions 0, § . 

Table 1: The ratio R = IZZ P /I™ 1 averaged over 30 nuclei. 



Models: 










R 


2.2(2) 


2.7(2) 


1.5(1) 


1.7(1) 



Precise 7-decay parameters are, however, necessary for the calculation of the interaction cross- 
sections of neutrons with unstudied target nuclei and the understanding of the behavior of nuclear 
matter with increasing excitation energy. An analysis of the existing methods for the determination 
of the level density || f?j and radiative strength functions (k) || 

k = r Xi /{E% x A 2 / 3 x D\) (1) 

in deformed nuclei, for example, shows that it is not possible to obtain sufficiently precise exper- 
imental level densities for certain intervals of their energies and quantum numbers as well as the 
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widths of the corresponding transitions. Analysis of contributions of different sources of system- 
atical errors in determination of the level density from the evaporation spectra was performed by 
H.Vonach || mainly for light spherical nuclei. The total uncertainty evaluated by him amounts to 
about 20-30%. It should be noted, however, that an accuracy in calculation of cascade intensities 
within the models ||, |[in the lightest nuclei (from the range of the minimum of the neutron 
strength function) studied by us is also considerably better than that for deformed nuclei: the ratio 
R for UA Cd and 124 > 125 Te varies in limits from 0.7 ( 124 Te) up to 1.4 ( U4 Cd). 

Therefore, without developing new methods for the determination of nuclear parameters under 
discussion one cannot expect any progress in the modification of the existing theoretical models, 
first of all, for deformed nuclei. (In eq. (|l|) T X i is the partial radiative width of 7-transition with 
the energy E 7 , D x is the average level spacing of the decaying state and quantity A is the nucleon 
number). 

A new and sufficiently perspective way to obtain such information for the entire energy interval 



below B n seems to be the investigation |lC], 11] of the two-step 7-cascades between the compound 
state A and the given low-lying level / through all possible intermediate states i. The algorithms 
12, 13, |l4j] developed for the analysis of 7 — 7 coincidences registered by ordinary Ge detectors 



allow one to determine the intensity distribution of the cascades as a function of the energy of 
the cascade intermediate levels over the whole energy region up to E ex ~ B n with an acceptable 
systematic error (which decreases as the efficiency of the 7-spectrometer increases). 



The intensity i 77 of an individual cascade is 

^ x ^, (2) 



where T X i and Tif are the partial widths of the transitions connecting the levels A — ► i — > /, T x 
and are the total widths of the decaying states A and i, respectively. The sum intensity I 77 of 
the cascades is related to an unknown number of intermediate levels n X i = p x AE and unknown 
widths of primary and secondary transitions via the equation 

j-f i T x Ti j^<T Xi >m Xi <T if >m if 

The summation is over a certain set of quantum numbers of intermediate, initial, and final 
states for the purpose of comparison with the experimental data. The thermal neutron capture 
cross-section for two possible spins of compound states are listed in p5| , for example. The J n values 
for the initial and final cascade levels are also known. The latter, however, is true if the energy Ef 
of the final state does not exceed ~ 1 MeV. The optimal width of the interval AE and the number 
N of such intervals in eq. (§) are determined by the statistics of 7 — 7 coincidences (as a square 
detector efficiency) and the necessity to obtain detailed energy dependence for J 77 . The width of 
AE does not exceed 0.5 MeV even in the case of a 10% efficiency detector, however. The total 
radiative widths T\ of the capturing states are also known from the corresponding experiments for 
all stable nuclei [15]. The mean partial widths < T Xi >, < Tif > and the total numbers m Xi , rriif 
of levels excited by El and Ml transitions after the decay of the states A and i, respectively, to be 
found in the analysis are related to the total radiative widths as 

T x = < T Xi > xm Xi 

Ti = < T if > xmif (4) 

The contribution of higher multipolarities to eqs. (0) and (Q) is smaller than the error of the 
determination of / 77 . Equations (||) and (||) and their obvious combination 

T A x J 77 = ^T Xi x rii x (T if / < T if > m if ) (5) 
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allow three ways of the estimation of the parameters of the cascade 7-decay using the experimental 
data on I 77 and T\: 

(a) the level density can be estimated from eq. (||) using model calculated partial radiative 
widths; 

(b) the partial widths of cascade transitions can be estimated from eq. (pj) using model 
calculated level densities with certain J 71 "; 

(c) simultaneous estimation of the intervals of probable level densities and radiative strength 
functions which satisfy eqs. ((^) and (Q) in general. 



It is clear that the level density and strength functions found according to variants (a) and (b) 
inevitably contain errors caused by the uncertainties of experimental and model values used as 
parameters of the analysis. However, the influence of these uncertainties on the final result is 
suppressed because of the correlation (determinated by the used type of the functional relations 
(|3|) and (||)) between the experimental T e ^ p , I^ p and the parameters under study p, T. 

In accordance with the variant (a) the sufficiently narrow interval of probable p was determined 
for almost 30 nuclei from the mass region 114 < A < 200 for some set of possible models of 7- 
transition strength functions. An important conclusion made in [16| is that the best description of 
the level density in the interval from ~ 0.5B n to B n was achieved in the framework of the general- 
ized model of the superfluid nucleus Q. Besides, simple enough models |B| of radiative strength 
functions cannot provide a correct description of the experiment and also need modification. An 
analysis by variant (b) was performed by us, as well. The main result is that there are no strength 
function models for El and Ml transitions in deformed nuclei which could reproduce the depen- 
dence T\ x J™ at primary transition energies E\ < 2 — 3 MeV if the level density is set by the 
model of a non- interacting Fermi-gas. Therefore, the understanding and correct description of the 
7-decay of the compound nucleus with a high level density require experimental determination of 
the level density and radiative strength functions over the entire excitation energy region. 

Further investigations [17| have shown that the level density at excitations from 1-2 to 3-4 MeV 
in, first of all, deformed nuclei deviates strongly from the exponential energy dependence derived 
on the basis of the idea that the nucleus is a non-interacting Fermi-gas [|J. Moreover, it is not 
excluded that the level density in this energy interval can be almost constant or even decrease with 
increasing excitation energy. These confirm and complement the results obtained in [16]. 



2 Analysis 

The variant (c) of analysis of the experimental intensities of two-step 7-cascades between the 
capturing state and several low-lying levels allowed us to suggest an original method for the solution 
(although partial) of this problem. It is based on an obvious circumstance that N + 1 equations 
@ and (||) together with 6A^ conditions 

p(n = +) > 0; p(ir = -) > 

T(E1) > 0; T(Afl) > (6) 

(separately for primary and secondary transitions in the case of radiative widths) restrict some 
interval of possible level densities and partial radiative widths which provide a simultaneous re- 
production of Y e ^ p and I 77 p . This interval can be estimated using modern computers and the 
existing computational algorithms. Its width, however, cannot equal zero even at zero uncertainty 
of the experiment. It should be added that / 77 in the form of eq. (|3|) is inversely proportional 
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(qualitatively) to the total number of states excited in the process under study and is proportional 
to the ratio of cascade transition widths to their mean values. Therefore, the method of analysis 
described below has a maximum sensitivity at minimum density of the excited states (unlike the 
methods §, §). 

As in the case of other reactions (followed by 7-emission) used for the determination of p, 
all values obtained experimentally in the (nth, 7) measurements are determined by the product 
Y\i x p. Hence, in the calculation deviation of one of the two parameters from its mean value 
is compensated by deviation of the other one with the corresponding magnitude and sign. This 
circumstance should be taken into account in data processing — a minimum or maximum value of 
the level density derived from the experimental data results, e.g., in a maximum or minimum value 
of the corresponding strength functions. 

It should be noted that deviation of the calculated level density from the true value is completely 
compensated by deviation of strength functions when T\ is only calculated. In the case of the 
calculation of J 77 the compensation is incomplete. This very circumstance allows one to select 
the intervals of p and T\i which provide the description of the / 77 and T\ parameters with an 
acceptable uncertainty. This analysis can be performed by means of finding large enough sets of 
random values of p and T\i which reproduce completely the parameters Y e ^ p and and belong 
to the intervals that contain true values. This means that most probable values of the level density 
and radiative strength functions of dipole 7-transitions and intervals of their uncertainties can be 
found by selection of pairs of random p and k which satisfy, in general, eqs. (||) and @ or (||) and 
(j5|). This requires numerous repetitions of the procedure and statistical methods of analysis. 

It is clear that the widths of the intervals of probable p and k satisfying eqs. (||) and (|j) in- 
crease with increasing number of unknown parameters in the equations. According to experimental 
conditions, the summation in eqs. (||) and (||) as over all intermediate states of the cascades. Since 
the summed data included cascade transitions of different multipolarities, we could not obtain the 
strength functions of El and Ml transitions and the level density for different parities separately 
with a good precision. In practice, from a combination of eqs. (0) and (|5|) the sum of strength 
functions and the sum of level densities of both parities should be only derived and compared 
with model predictions. The corresponding summation reduces considerably the uncertainty of the 
observed result due to anti-correlation of elements. 

Indeed, an analysis of the available data confirms that the dispersion of each set of p(ir = +), 
p(-K = — ), k(El) and k(Ml) random values is too large to make any conclusions about independent 
correspondence of individual values to the model. 

A sufficiently large N and the nonlinearity of eqs. (||) and (|j) stipulate the choice of the way 
to solve the system of equations and inequalities - the Monte Carlo method. The simplest iterative 
algorithm [jl6| was used for this aim: we set some initial values for T(E1), T(M1), p(ir = —), and 
p(ir = +) and then distort them by means of random functions. If these distortions decrease the 
parameters A = (J 77 p — I 77 ) 2 at this step of the iteration procedure, then the distorted values 
are used as initial parameters for the next iteration. Agreement between the experimental and 
calculated cascade intensities and the total radiative widths, respectively, is usually achieved after 
several thousand iterations. As a result we get two random ensembles of level densities and partial 
widths for every N energy intervals. Examples of intermediate and final results of one of many 
variants of the calculation for two nuclei are shown in Figs. 1 and 2. It is obvious that such iterative 
process can be realized in an unlimited number of ways. We chose a sufficiently simple and effective 
way: the Gaussian curve is used as a distorting function for logarithms of p and / 

f(E) = Axexp(-(E-E ) 2 /a 2 ) (7) 

Its parameters are independently chosen for the level density and strength functions from the 
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intervals [-0.2;0.2], [E^; B n ] and [0.3 MeV;B n ] for A, E , and a, respectively using a standardized 
random value distributed uniformly in [0;1]. Here Ed is the maximum excitation energy of the 
known discrete level involved in the calculation. Numerous repetitions of the iterative calculation 
with different initial parameters (including obviously unreal values of T and p) for ~30 nuclei from 
the mass region 114 < A < 200 show that this algorithm yields rather narrow intervals of the sum 
level density of both parities and of the sum partial widths of El and Ml transitions. The use of 
eq. (0) with mentioned parameters allows one to get a set of different, smooth enough functional 
dependences for both p and k. In this case for the majority of the studied nuclei the values of level 
density are in good agreement with the number of the observed intermediate levels of the cascades 
resolved as the pairs of peaks. In some nuclei, however, the mean level density (which together with 
the mean strength functions provides reproduction of cascade intensities) is less than the number 
of intermediate levels observed below ~ 2 MeV. The main portion of this discrepancy is removed in 
all cases if one foresees a possibility of additional local variation of k for high-energy transitions in 
the energy interval which, as a rule, does not exceed 0.1-0.2 MeV. One of the examples of this kind 
is shown in Fig. 1. The necessity to account this effect can be due to both insufficient averaging of 
the random partial widths of primary transitions and their possible dependence on the structure 
of the excited low-lying level. This can result, for instance, from concentration of the strength of 
the fragmented single-particle or phonon states. 

3 Asymptotical uncertainty of the obtained parameters 

The method suggested by us for estimation of p k cannot give unique value of these parameters 
at a given energy of the excitation or quantum energy. Therefore the question arises about the 
value of their uncertainty at different energies and degree of possible systematical deviations of 
the observed parameters from the modal values. The results of modelling for 15& Gd and 198 ^4u 
shown in Fig. 3 answer these questions. Intensity of cascades for these nuclei were calculated under 
assumption that the strength function k(El) is described by model Q and value of k(Ml) = const; 
level density exponentially increases with the energy or have some step-like structure. Below the 
excitation energy ~ 1 — 2 MeV the calculation used experimental decay scheme. Consequently, the 
calculated intensity distribution of cascades in function of the primary transition energy has one 
or two maxima. (Other conditions of the calculation completely corresponded to the experiment). 

Figure 3(b) shows that the model level density is reproduced practically without systematical 
error and the width of the interval of its probable values does not exceed 20-30%. 

Discrepancy between the experimental and model sum k{El) + k(Ml) results from that the 
total radiative width calculated according model 0] does not correspond to the experimental value. 
Energy dependence of k(El) + k(Ml) is reproduce rather well - sharp changes in the first derivative 
with respect to the quantum energy is not observed (unlike some other nuclei studied by us). 

So, one can summarize that the suggested method provides reliable enough estimation of the 
level density and radiative strength functions of dipole transitions. 

4 Approach used in calculation 

The insufficient experimental data on cascade 7-transitions (only cascades terminating at low-lying 
levels (Ej < 1 MeV) of nuclei were studied H) does not allow us to determine the level densities and 
gamma-widths without the following important assumption: the strength functions of transitions 
of a given multipolarity only depend on the transition energy and do not depend on the structure 
and energy of the corresponding excited states. Their nonequal values for 7-transitions of equal 
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energies but populating different levels is, in part, compensated by the circumstance that the left 
part of eq. (|5]) depends on absolute radiative strength function values of primary transitions and 
depend only on the ratio of strength functions in the case of secondary transitions. These decrease 
the effect of the discussed assumption on the k(El) + k(Ml) values but do not remove it completely. 
There is no necessity in introduction of any hypotheses of spin dependence of level density differing 
from that predicted in models |^, ||] . 



5 Sources of errors in the determination of strength functions 
and level densities 

The presence of the statistic and systematic errors in determination of J~~, T\ and specific problems 
of extraction of level density and radiative strength functions cause noticeable uncertainties of the 
determined parameters. The influence of the different sources of errors on the obtained results 
manifest itself in a different degree. 

1. Uncertainties of the measuring of terms in eqs. (|3|) and (||) result in errors of strength functions 
and level density. Owing to a linear relation between T\, I 77 and F\i in eq. (||), ~ 10% errors 
of T\ and i™ achieved in the experiment cause rather a small error in the determination of 
I\j and p as compared to dispersion of the obtained data. 

2. The more considerable source of uncertainty in the determination of the strength functions 
and p is a systematic error of decomposition [jD|, [l4j] of the experimental spectra into two 
components corresponding to solely primary and solely secondary transitions. The analysis 



| I3I showed that the error in AI 77 caused by this procedure does not usually exceed ~ 20% for 
primary transition energy E\ < 3 — 4 MeV. Intensities of cascades (histograms in Figs. 1, 2, 
4-13) at these primary transition energies can be overestimated, as a maximum, by the above 
value. At the higher energies they can be decreased by the same value (the total intensity 
is preserved). In order to estimate the influence of A/ 77 on the final results, the I 77 values 
were varied within a level of 25%. These variations caused changes in k{E\) + k{M\) and p 
which did not exceed the dispersion of the data plotted in Figs. 4-13. 

3. The maximum uncertainty of level density and radiative strength functions results from the 
use of condition (6). It dominates at any possible precision in determination of i" 77 and T\. 
The simplest way to estimate these errors at any E\ and E ex is the following: 

(a) taking into account that the probabilities of deviations with opposite sign of the random 
Pi and fi values with respect to their mean values are equal and decrease as the absolute 
values of deviations increase; and 

(b) assuming that mathematical expectations of the random ensembles of the pi and ki 
values satisfying eqs. ©-(H) correspond to their real values 

one can consider the mean-square deviations of the random values relative to their arithmetical 
means as the estimations of the errors. These errors can be attributed to level density and 
strength functions separately in spite of their strong anti-correlation. Just these uncertainties 
are shown for the radiative strength functions and level density plotted in Figs. 4-13. 

On the whole we can summarize the situation as the following. At the presently achieved 
accuracy for experimental determination of i" 77 and T\, level densities and strength functions 
are derived from eqs. with the mean total uncertainties of about 40-50% in the worst 

case. Asymptotic value of this uncertainty at zero statistic and systematic errors of the 
experiment is equal, in the average, for both p and k and cannot be less than ~ 20%. 
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4. There are two ways to decrease the errors of the level density and strength functions deter- 
mined from eqs. (H)-©: 

(a) the increase of the volume of the experimental data on the cascade intensities; 

(b) the reduction of the number of parameters in eqs. owing to the use of additional 
information or introduction of some new assumptions. 

In the first case the problem can be easily solved experimentally: the use of a Compton- 
suppressing spectrometer consisting of HPGe detectors with an efficiency of not less than 
30-40% allows the selection from a mass of 7 — 7 coincidences of two-step cascades for a 
considerably larger number of their final levels than at present. From a combination of eq. (^) 
for the sum over all final levels of cascades and an individual final level / one can determine 
the ratio Y^f/ < T^ > 'xm i f for all possible values of i and /, i. e., determine energy 
dependence of the experimental sum k{E\) + fc(Ml) for any possible secondary transitions, 
get rid of the only approach used in the analysis, and reduce the number of parameters in 
the analysis. 

The data shown in Figs. 4-13 were obtained under assumption about a constancy of the ratio 

Y xi/Y if = const (8) 

for the transitions with equal multipolarity and energy E\ in all interval of the neutron binding 
energy. 

The comparison of the total 7-spectra and population of low-lying levels calculated in this way 
with the available experimental data including spectroscopic information |jj shows that even such 
assumption provides better accuracy in calculating the parameters of cascade 7-decay than the 
approach using the models [2|, [3], ||, ||. Unfortunately, we could not achieve complete correspon- 
dence between the estimated level density and available spectroscopic information. Nevertheless, 
the obtained values demonstrate certain correspondence of our level density with the numbers 
of the excited levels observed in the experiment jlj. Some nuclei, however, demonstrate residual 
discrepancy (for example, 170 Tm, Fig. 1). This discrepancy can be attributed, partially, to both 
insufficient precision of assumption (|8|) and inexactitude of the spectroscopic data. Their errors can 
be also considerably decreased using more efficient spectrometer of 7 — 7 coincidences than that 
used by authors 0. 

On the whole, in spite of the uncertainties mentioned above one can conclude that at a given 
stage of the experimental investigation of the cascade 7-decay of compound states our method 
provides more reliable results than methods || [/], ||] . 

6 Main results of analysis 

The type of relation between k and p on the one hand and between Y\ and I 77 on the other hand 
does not allow one to determine k and p unambiguously and independently. Some deviation of, for 
example, p from a real value is inevitably compensated by deviation of strength functions of the 
corresponding magnitude and sign. Nevertheless, the results obtained in the present analysis can be 
used for the verification of nuclear models and, if necessary, for the determination of the direction 
of the further development of these models. The main argument in favour of this statement is 
relatively week dependence of the final results on the initial values of strength functions and p 
in the iterative process. As an example, Figs. 1 and 2 show the strength function and p values 
obtained for their unreal initial values: p{E ex ) = p{B n ), the strength functions decrease linearly as 
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the transition energy increases. Nevertheless, the final results of the iterative process quite agree 
with a general picture obtained for a large enough set of different real and unreal initial values of 
k and p. This confirmes the conclusion that the strength functions and level density obtained from 
the analysis can be considered as most probable. 

The strength functions k(El) + k(Ml) and level densities p obtained in the present analysis are 
plotted in Figs. 4-13. For every set of random p at a given excitation energy E ex and k(El) + k(Ml) 
at a given primary transition energy E\ = B n — E ex there were determined both their mean values 
and probable dispersion using usual relationships of statistical mathematics. The results of the 
analysis are compared with predictions of the level density models (2|, ||] and models of radiative 
widths [§, H . In the case of radiative strength functions a comparison is performed in the following 
manner: the k(El) values calculated according to the models j| and ]|] (upper and lower curves, 
respectively) are summed with k(Ml) = const which is normalized so that the ratio T(M1)/T(E1) 
would be approximately equal to the experimental data at E^ ~ B n . 

A comparison of the results of the analysis with predictions of the models ||, |3|, ||, || (often used 
by experimentalists) shows that: 



1. the sums k{El)+k{M\) and p are not monotonic functions of the energy and, probably, reflect 
the most common peculiarities of the structures of the states connected by the corresponding 
7-transitions; 

2. the energy dependence of k(El) + k(Ml) differs strongly from predictions of the models [|], |5| 
in the case of even-even compound nuclei from the region of the 4s-resonance of the neutron 
strength function, at least; 

3. the k(El) + fc(Ml) functions increase from near-magic to deformed nuclei and from compli- 
cated highly-excited states to simpler low-lying levels which are populated by 7-transitions 
under consideration; 

4. relative deviations of the obtained strength functions and level densities from the mean values 
are characterized by strong negative correlation. In the majority of nuclei the correlation 
coefficient changes from -0.6 to -1.0. This means that the strength functions and level densities 
are not independent variables in eqs. (||) and (|5|), which provides the possibility of their 
simultaneous determination; 

5. the probable level density determined in the present analysis conforms to the picture obtained 



in previous experiments [16, 17]: up to the excitation energy 1-2 MeV, our data are not in 
contradiction with the exponential extrapolation of p{E ex ) predicted by the Fermi-gas back- 
shift model 0. The energy dependence of the level density in the interval from 1-2 to some 
threshold value £"& is considerably weaker than it follows from any existing level density 
model. Above Ef, ~ 3 MeV for iV-odd and ~ 4 MeV for iV-even nuclei, the level density, 
most probably, corresponds better to the predictions of the generalized model of the superfluid 
nucleus in its simplest form |J. 



This change in the behaviour of the level density in the vicinity of the excitation energy E\, may 
signify a qualitative change in the nuclear properties. The observation [18] of the probable har- 
monicity of the excitation spectra of the intermediate levels of the most intense cascades in a large 
group of nuclei from the mass region 114 < A < 200 allows an assumption that the nuclear prop- 
erties at low energy are mainly determined by vibrational excitations (probably, a few phonons of 
rather high energy). A very quick exponential increase in the level density above E^ says about the 
probable dominant influence of the inner, many-quasiparticle type of excitations of these states. 
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7 Discussion 



The method suggested in present work allows model independent, simultaneous estimation of in- 
tervals of probable values of the level densities with given spins and summed strength functions of 
primary dipole transitions populating them. The method is effective in investigations of any stable 
nucleus. The main differences of this algorithm from the known methods of determination of level 
densities || [7| [| and radiative strength functions [|| are the following: 



1. Our method does not permit one to get the sole values of p and k for a given energy. But the 
width of the intervals of their probable magnitudes depends very weakly on the uncertainty 
in determination of T\ and I~~ at the achieved precision of the experiment, at the one hand, 
and is narrow enough in order to get new information on nuclear matter, from the other hand. 

2. The most correct and reliable data on the level density is derived from the evaporation spectra 
at the highest excitation energies; analysis of the cascade intensities provides similar data for 
the lowest energies. So, both methods mutually add each other. 

3. Analysis of cascade intensities allows direct determination of the absolute level densities, 
evaporation spectra usually provide || f?J information on nuclear temperature. 

4. Systematical uncertainties of both methods do not relate. Discrepancies in the independently 
determined level densities at some energies indicate to necessity, for example, to determine 
more precisely the barrier transmission factor for the evaporated particle or to take into 
account different energy dependence of k of the primary and secondary transitions of the 
7-cascades. Besides, they can testify to necessity to describe more correctly direct and pre- 
equilibrium processes in nuclear reactions for deformed nuclei or to define more precisely 
the nuclear excitation energy above which thermodynamical parameters of a nucleus are 
determined mainly by quasiparticle excitations. 

5. Energy dependence of the data in Figs. 4-13 can be reproduced well enough in the framework 
of modern version of the generalized model of the superfluid nucleus ]l9| if the temperature 
of the phase transition is diminished up to the value T' cr ~ 0.7T cr , where 

T cr = 5/1.76 (9) 

is the temperature of the transition from the superfluid to normal phase of homogeneous 
Fermi-system pQ]. But re-determination of the entropy and temperature predicted by model 



19] should be done so that nuclear temperature below T' cr will not increase with decreasing 



excitation energy. 

6. Additional and independent arguments in favor of reliability of step-like structure in level 
density are: 

(a) combinatorical calculation |^TJ of density of the states with K n = 1/2+ in 1G5 Dy below 
B n , providing similar to Figs. 4-13 picture; 

(b) analysis (22] of the experimental data from the reaction 1&5 Ho(p,n) ie5 Er. This also 
demonstrates some step-like structure in the total level density at low excitations; 



(c) precise analysis [23| of the neutron cross sections for actinides testify to necessity to take 
into account the influence of the pairing interaction on the level density for the wide 
interval of the neutron energies manifesting itself, in particular, as irregularities in the 
energy dependence of the level density. 
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7. It is obvious that the structures shown in Figs. 4-13 can be inherent not to the total level 
density with given J 71 ", but only to that part of them which are really excited in (n, 7) reaction. 
Then, unlike the existing notions, this reaction is selective and structures of the excited states 
must be taken into account in any calculations of parameters of this reaction in the entire 
excitation energy region below B n . 



8 Conclusions 

A new method is suggested for a simultaneous estimation of the probable level density populated by 
dipole primary transitions in the (nth, 7) reaction and the sum strength functions k(El) + k(Ml) of 
these transitions. Unlike other methods used for the investigations of nuclear properties below the 
excitation energy 6-9 MeV, this method allows the estimation of p, radiative strength functions, 
and intervals of their probable variations without any model notions of the nucleus. 

The method is universal - it can be used for any nucleus and reaction with 7-emission. The 
latter is possible if the excitation energy interval of high-lying states is narrow enough in order to 
use the sum coincidence technique. Besides, the most probable quanta ordering in the cascades 
must be determined for the main part of the observed cascade intensity. It should be noted, that 
in the case of a lack of the experimental values of the total radiative widths of decaying high-lying 
states the absolute radiative strength functions cannot be determined. In this case only relative 
energy dependence of the radiative strength functions can be obtained. 

The most important (although preliminary and qualitative) physical result is that the level 
density below the neutron binding energy (first of all in deformed nuclei) cannot be reproduced 
to a precision achieved in the experiment without more precise than in accounting for the co- 
existence and interaction of superfluid and usual phases of nuclear matter in this whole excitation 
energy interval. 

The obtained results demonstrate very serious and obvious discrepancies with the existing 
ideas of the structure of the deformed nuclei. These data agree completely with an earlier obtained 
qualitative picture IS] of the studied process: considerable influence of vibrational excitations on 



the nuclear properties below the excitation energy Eb and a transition to dominant influence of 
quasiparticle excitations above this energy. 

This work was supported by RFBR Grant No. 99-02-17863 



10 



References 

[I] S.T. Boneva et al., Phys. of Atomic Nuclei 62(5) (1999) 832 

[2] W. Dilg, W. Schantl, H. Vonach, M. Uhl, Nucl. Phys A217 (1973) 269 

[3] A.V. Ignatyuk, Proc. of IAEA Consultants Meeting on the use of Nuclear Theory and Neutron 
Nuclear Data Evaluation (Trieste, 1975) IAEA-190, Vol. 1 (1976) P.211 

[4] P. Axel, Phys. Rev. 126(2) (1962) 671 

J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear Physics, New York (1952) 

[5] S.G. Kadmenskij, V.P. Markushev, W.I. Furman, Sov. J. Nucl. Phys. 37 (1983) 165 

[6] M.I. Svirin, G.N. Smirenkin, Yad. Fiz. 47 (1988) 84 

[7] A.V. Ignatyuk et al., Sov. J. Nucl. Phys. 29 (1978) 875 

[8] G.A. Bartholomew et al., Advances in nuclear physics 7 (1973) 229 

[9] H. Vonach, Proc. IAEA Advisory Group Meeting on Basic and Applied Problems of Nuclear 
Level Densities (New York, 1983), INDC(USA)-092/L, (1983) P.247 

[10] S.T. Boneva et al., Sov. J. Part. Nucl. 22(2) (1991) 232 

[II] S.T. Boneva et al., Sov. J. Part. Nucl. 22(6) (1991) 698 

[12] Yu.P. Popov, A.M. Sukhovoj, V.A. Khitrov, Yu.S. Yazvitsky, Izv. AN SSSR, Ser. Fiz. 48 
(1984) 1830 

[13] S.T. Boneva, V.A. Khitrov, A.M. Sukhovoj, A.V. Vojnov, Z. Phys. A338 (1991) 319 

[14] S.T. Boneva, V.A. Khitrov, A.M. Sukhovoj, A.V. Vojnov, Nucl. Phys. A589 (1995) 293 

[15] S.F. Mughabghab, Neutron Cross Sections V.l, part B, NY.: Academic Press (1984) 406 

[16] V.A. Khitrov, A.M. Sukhovoj, Proc. of VI International Seminar on Interaction of Neutrons 
with Nuclei (Dubna, 1998) E3-98-202, Dubna, 1988, P. 172 

[17] A.M. Sukhovoj, V.A. Khitrov, Physics of Atomic Nuclei 62(1) (1999) 19 

[18] A.M. Sukhovoj, V.A. Khitrov, Bulletin of the Russian Academy of Sciences, Physics 61(11) 
(1997) 1611 

[19] A.V. Ignatyuk, Report IAEA INDC-233(L) IAEA, Vienna 1985 

E. M. Rastopchin, M. I. Svirin and G. N. Smirenkin, Yad. Fiz 52, (1990) 1258 (in Russian) 

[20] J. Bardin, L. Cooper, J. Schriffer, Phys.Rev. 108 (1957) 1175 

[21] A.I. Vdovin et al., Sov. J. Part. Nucl. 7(4) (1976) 380 

[22] B.V. Zhuravlev, Izv. RAN, Ser. Fiz. 63(1) (1999) 148 (in Russian) 

[23] V.M. Maslov, Phys. of Atomic Nuclei 63(2) (2000) 161 



11 




Fig. 1. The examples of p (a) and k (b) intermediate values and the corresponding distributions of 
cascade intensities (c) for the 170 Tm odd-odd nucleus in function of the primary transition energy 
E\ or excitation energy E ex . Letters next to the lines mean the number of iterations. Triangles 
show number of levels excited by the primary dipole transitions with the energy E\ in the energy 
interval of 100 keV. The dashed curve (a), (b) represents model predictions, the histograms (c) 
represent the experimental cascade intensities with statistical errors. 
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2. The same as in Fig. 1, for the 190 Os even-even nucleus. 
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Fig. 3. The intensity of two-step cascades (a) calculated using level density [| shown by solid lines 
in (b) and radiative strength [Q, ||] functions - line 1 in (c) (line 2 in (c) represents predictions 
of model [5]). Points with error bars represent the interval of possible values of p (b) and k (c) 
providing acceptable precision in reproduction of cascade intensities shown in (a). 
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Fig. 4. Experimental cascade intensities I~ n in 0.5 MeV energy bins with ordinary statistical errors 
for 114 Cd and Te (histograms. Curves represent calculation performed like that shown in Fig. 3. 
Points with errors represent number of levels per 100 keV energy interval and sums k(El) + k(Ml), 
respectively. 
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Fig. 5. The same as in Fig. 4 for 137 Ba and 138 Ba. 
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Fig. 6. The same as in Fig. 4 for 139 Ba and 146 Nd. 
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Fig. 7. The same as in Fig. 4 for 150 Sm and 156 Gd. 
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Fig. 8. The same as in Fig. 4 for 15S Gd and 160 T6. 
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Fig. 9. The same as in Fig. 4 for 164 Dy and 168 Er. 
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Fig. 10. The same as in Fig. 4 for 170 Tm and lu Yb. 
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11. The same as in Fig. 4 for 181 Hf and 183 W. 
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Fig. 12. The same as in Fig. 4 for 192 Ir and 196 Pt. 



23 




Fig. 13. The same as in Fig. 4 for 198 Au and 200 Hg. 



24 



